Perfect imaging with positive refraction in three dimensions 



ON 

o 
o 

(N 
> 

o 



o 

H— > ■ 

o : 

o 

Ph. 



> 
(N 

in 
in 
o 



ON 
O 



X 



Ulf Leonhardt and Thomas G. Philbin 
University of St Andrews, School of Physics and Astronomy, North Haugh, St Andrews KY16 9SS, UK 

(Dated: November 3, 2009) 

Maxwell's fish eye has been known to be a perfect lens within the validity range of ray optics since 
1854. Solving Maxwell's equations we show that the fish-eye lens in three dimensions has unlimited 
resolution for electromagnetic waves. 
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Perfect imaging with positive-index materials has been 
discussed since 1854 but only for light rays and not 
for waves @, Q. In fact, the defining property of such 
ideal optical instruments Q is the perfect focussing of 
rays: for each point in a region of space, all emitted rays 
meet in a corresponding point in the image space. Yet 
it is the wave nature of light that sets the resolution of 
optical instruments. Only very recently one of us 
proved that the archetype of all ideal optical instruments 
[2| , Maxwell's fish eye [Ij perfectly images light waves in 
two dimensions. Here we show that the fish eye in three 
dimensions has infinite resolution as well, provided its 
medium is impedance-matched. Absorption does not ap- 
pear to significantly reduce the image quality, in contrast 
to imaging with negative refraction [5j, |6j. However, fish- 
eye lenses contain both the source and the image inside 
the medium and impedance-matched devices are still dif- 
ficult to make in practice. On the other hand, our case 
supports the idea |4| that perfect imaging is not neces- 
sarily caused by the amplification of evanescent waves [f| 
but rather by the geometry of light 0, H, [1] • Such con- 
ceptual insights may be vital for the future direction of 
technical developments for perfect imaging. 

According to Fermat's principle @, Q, light rays fol- 
low extremal optical paths with a path length measured 
by the refractive index n — geodesies. Optical mate- 
rials thus establish virtual spatial geometries for rays 
@, Ho} - Whether these geometries are also valid for 
waves depends on the type of wave propagation. In two- 
dimensional structures, typical of integrated optics, the 
polarization of light decides whether the structures are 
perceived as geometries: for materials with purely electri- 
cal response TE polarization is required [i| . In three 
dimensions, the material must be impedance-matched 
for establishing a virtual geometry for electromagnetic 
fields 0] , with equal electric permittivity e and magnetic 
permeability /z, 



(1) 



This crucial condition has not been considered in the 
previous treatment of Maxwell's fish eye with Maxwell's 
equations [H, EH- Scalar waves [H should not obey 
the Helmholtz or Schrodinger equation either, for perfect 
imaging in three dimensions, but rather a wave equation 
we consider here as well, Eq. (fTpl) In t he p revious wave 
theories of Maxwell's fish eye [l3, El El El, perfect 



imaging was impossible. If, however, the medium or po- 
larization is chosen such that the geometry of light is not 
restricted to rays but extends to waves, waves may be as 
perfectly imaged as rays 0]. 

Maxwell's fish eye turns out 0] to establish a non- 
Euclidean geometry [lj|, the 3D surface of the 4D hy- 
persphere, by the refractive-index profile 



(2) 
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where r denotes the distance from the center measured in 
the characteristic length scale of the device. As the index 
profile ([2]) is radially symmetric, the trajectory of a light 
ray lies in a plane, due to the conservation of angular 
momentum. So, for light rays, the propagation in 3D is 
the same as in 2D where Maxwell's fish eye corresponds to 
the surface of a 3D sphere. Here all the rays emitted from 
one point travel along the great circles, meeting again at 
the antipodal point, which proves that Maxwell's fish eye 
fits the definition of an ideal optical instrument @. In 
physical space, the antipodal image of a point r turns 
out [3J] to appear at 



ro 

" ^2 



(3) 



In this paper we show that electromagnetic waves are 
perfectly imaged at r , apart from a phase delay of 



if = nk 



(4) 



where k denotes the wavenumbcr in units of the inverse 
length scale of the fish eye. As the phase (0| is uniform, 
objects are not only faithfully but also coherently imaged. 
Furthermore, the phase ^ is linear in wavenumber and 
so in frequency if the refractive index is not frequency- 
dependent. In this case, the time delay between source 
and image is uniform as well. Note that the phase de- 
lay is different in 2D [4|: ttv with v 2 = k(k + 1). In 
the following we prove perfect imaging and obtain the 
quantitative results ([3]) and (|4]) by analyzing the electro- 
magnetic Green's function. 

Electric and magnetic fields. — The Green's function 
G describes the electric field of a stationary wave with 
wavenumber k emitted by an elementary dipole at po- 
sition ro that may point in all three directions. The 
Green's function is a matrix also known, in electrical en- 
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gineering, as a dyade or a bi-tensor with the first in- 
dex referring to the electric-field strength at the spec- 
tator point r and the second describing the direction of 
the dipole source at ro. If not otherwise stated we use 
Cartesian coordinates. As any source is a weighted and 
directed distribution of elementary point dipoles, it is 
sufficient to consider the Green's function G(r,ro) for 
perfect imaging. The Green's function obeys the wave 
equation 



V x -V x G 

n 



nk 2 G = it 



(5) 



with the infinitely localized current density i. In a sta- 
tionary regime, the wave emitted at the source must dis- 
appear somewhere. Usually the wave ultimately disap- 
pears at oo, but in perfect imaging the entire energy of 
the emitted wave is focused at the image and does not 
reach oo at all. In practice, the light will be absorbed at 
the image, for example in a lithographic photoresist or a 
detector. In any case, in order to maintain a stationary 
regime we must supplement the source by a drain at the 
image Q with phase delay (Tj|J) such that [4| 



6(r - r ) - e^5(r - r' ) . 



(6) 



The minus indicates that the drain is a source run in 
reverse. Finally, causality implies [4] that G is analytic 
in k on the upper half complex plane and vanishing for 
k — » oo there. If G meets all these requirements our case 
is settled: Maxwell's fish eye makes a perfect lens in 3D. 

In order to deduce the Green's function it is advanta- 
geous to represent G in terms of its magnetic field, the 
tensor H , as 



G = 



V x H -it 

n k 2 



where we require H to obey the wave equation 

1 , it 

V x -V x H -nk 2 H = V x — . 
n n 



(7) 



(8) 



From Eqs. {7} and (J5J) follows Faraday's law of induction 



V xG = nH 



(9) 



and subsequently the defining wave equation (O of the 
Green's function. Faraday's law also reveals that H 
describes the magnetic field strength divided by ic/xofc 
where c denotes the speed of light and no the permeabil- 
ity of vacuum. Consider a special case of wave propaga- 
tion first. 

Source at origin. — Imagine the source is placed at the 
center of the fish eye (and the image would be at oo), 



ro = 0, 



6(r). 



(10) 



Similar to the electromagnetic wave emitted by a point 
dipole in free space [l7|, we represent H by the ansatz 



H = V x 2D(r) 1 



(11) 



In free space, 2D describes the scalar Green's function 
[17]. In Maxwell's fish eye, D turns out to be the scalar 
Green's function, as we will show, where the factor of 
2 represents the refractive index © at the source. We 
apply the expression 



V/(r) = - d r f 



(12) 



for the gradient of any radial function where d r abbrevi- 
ates the derivative with respect to the radius, and obtain 



H = r x 1 . 



(13) 



We calculate the curl of H by expanding the double vec- 
tor product in V x (V x 2Dt), using the radial gradient 
(fT2")l and the standard expression of the Euclidean Lapla- 
cian in spherical coordinates, 

V x H = V <g> ~ d r 2D -l(d?+ ^d^j 2D 

= ^(#-;fc)22>-l + 2D. (14) 

Then we turn to V x n~ 1 V x H that occurs in the wave 
equation {SJ of the magnetic field. For the first term 
in expression (|14p the only non-zero contribution to the 
curl originates from V x a ® r = —a x V (g> r = —a x 1. 
For the curl of the second term divided by n we apply 
the radial gradient (fT2|) as in formula (fl~3|) . Combining 
these terms we obtain 

1 2r / 2 1 
Vx-VxiJ= — xlf d r -d r r | d r D 



n 
2nr 

r 



X td r 



d r r 2 nd r D - D 



(15) 



where in the last step we used the explicit formula (|2|) 
of the fish eye profile. All expressions in the wave equa- 
tion (|5J), including the curl of the delta functions on the 
right-hand side, have the same matrix structure as the 
magnetic field (fl~3|) . We only need to require 

1 



r 2 71 3 



d r r 2 nd r D + (k 2 - 1)D = (16) 

4 n 

for finding the Green's function G in the special case 
(jTCTj) . But before we write down the solution of the radi- 
ally symmetric scalar wave equation (|16p we cast it in a 
geometric form and transform it, for investigating perfect 
imaging of scalar waves in 3D. 

Scalar waves. — The virtual geometry of the "fish-eye 
world" is characterized by the line element 



ds = n(r)dl . 



(17) 



In this geometry the Laplacian 7} appears in spherical 
coordinates r, 9 and <fi as 

E V - VQ = -^ d r r2n 9 r + \d 2 + ^^-d 2 



r*n° 
1 



rV • nV 



rr' 



(18) 
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in Cartesian coordinates. Consequently, we can write the 
radial wave equation (|16[) as 



-^V • nV 

n 6 



k 2 - 1 \D = =■ . 



(19) 



This is the wave equation of the scalar Green's function. 
The equation is invariant under coordinate transforma- 
tions that preserve the geometry, the line element (flTf . 
because it is entirely composed of geometrical construc- 
tions — the left-hand side contains the Laplacian and 
k 2 — 1, a scalar, while the right-hand side consists of 
delta functions divided by the density n 3 , the measure of 
length (fTT|) cubed. Such coordinate transformations are 
the Mobius transformations 



r(l 



r (l + 2r-r 



1 + 2 r • ro + r 2 r t 



2 r 2 



because one verifies that 



ds' = ds . 



(20) 



(21) 



In 2D, the Mobius transformations correspond to the ro- 
tations on the surface of the 3D sphere in stereographic 
projection on the plane 18]. In 3D, the Mobius transfor- 
mations ([2H]) describe the rotations on the 3D surface of 
the 4D hypersphere, the virtual geometry of Maxwell's 
fish eye [2|. As the relationship between source and 
image is also invariant under Mobius transformations, all 
scalar Green's functions are simply Mobius-transformed 
solutions of the radial wave equation (fTtJ)) with r replaced 
by the radius 



r-r„ 



v /l + 2r-r +r 2 7-2 



(22) 



A suitable solution is 
1 



D 



.r'-\ ) cxp(2ifc arctanr') 

8tt \ r 1 



If, 1 \ /1 + ir' 



8tt 



1 - ir' 



(23) 



a curious variation of the free space scalar Green's func- 
tion [l3]. The Green's function ([23]) has two singularities, 
one at r' — and the other at r' = oo; one singularity 
at the source point ro where the expression (|22|) van- 
ishes and the other singularity at the image (J3j) where the 
transformed radius ([22]) diverges. The phase 2k arctanr' 
is finite, because the "fish-eye world" is closed — the sur- 
face at the hypersphere is finite. At the image we obtain 
the phase delay j4}. The prefactor of the scalar Green's 
function ([23]) was chosen to give delta functions on the 
right-hand side of the wave equation (]T1?]) . The Green's 
function (|23[) is analytic in k and decays exponentially on 
the upper half complex k plane: D is causal. All our re- 
quirements are met: scalar waves with the wave equation 
([HI]) and the fish-eye profile @ support perfect imaging. 



Green tensor. — Having established the scalar Green's 
function D we follow a similar procedure to deduce the 
electromagnetic Green tensor G. First, we write down 
the Green's function for the source at the origin in 
Mobius-primed coordinates. One verifies that formula 
(0 with expression (fT4]) appears as 



G = 



V x n(r')V(g) V D(r')x V il 
n(r) n(ro) k 2 n k 2 



(24) 



evaluated at ro = 0. Here Vo denotes the gradient op- 
erator with respect to the source point ro and the arrow 
indicates that all terms on the left of it are to be dif- 
ferentiated. Second, all we need to do for establishing 
the Green's function for an arbitrary source point ro is 
Mobius-boosting the special case (fMj) . We know that r' 
is Mobius-invariant, but we also need to transform the 
bi-tensor components of G. For this, we write our result 
(]2~4"]) in geometric terms for the "fish-eye world" with line 
element (|17p . which is most easily done in index notation. 
We use the permutation symbol [abc] for the Cartesian 
curls [7[ and express them by two Levi-Civita tensors 
e abc g m fi sn _ e y e geometry, for spectator and source 
points separately. Finally we lower the first index of each 
Levi-Civita tensor by the corresponding metric tensors 
n(r) 2 l and n(ro) 2 l, respectively. In formulae, 



G a b — 



y [acd][bef] d 2 n(r') d 2 D(r>) iS ab 
n{r)n(r )k 2 dx c dxf, dx d dx{ nk 2 



cdef 



cd e ef 9 2 n ( r /) Q2 D ^ 

dx c dxi dx d dx f 



k 2 



l<5ab 

n k 2 



(25) 



Formula (|25|) describes a perfect bi-tensor in the fish-eye 
geometry that remains invariant after Mobius transfor- 
mation. Consequently, we can simply drop the quali- 
fication ro = 0: our result (|24p is valid for arbitrary 
source points. Note that the delta currents ((6]) in expres- 
sion (]2~4"]) conveniently cancel the delta functions arising 
in the derivatives of D; the Green's function G is sin- 
gular at source ro and image r' but does not develop 
delta peaks. As the scalar Green's function D is analytic 
and exponentially decaying on the upper half complex k 
plane, so is the Green tensor G. Source, image ^ and 
phase delay (H)) are the same as in the scalar case as well. 
In short, Maxwell's fish eye in three dimensions perfectly 
images electromagnetic waves. 

Mirror. — In practice, the fish-eye profile poses 
a formidable challenge: it is infinitely extended across 
space and the refractive index n < 1 for r > 1, the speed 
of light exceeds c here and approaches infinity. Both 
problems are related: Maxwell's fish eye represents a fi- 
nite virtual space, the 3D surface of the 4D hypersphere, 
stretched out to infinite physical space. Light can only 
reach infinity in a finite time with infinite speed. But, one 
can solve both problems in one stroke by placing a mirror 
at the unit sphere (r = 1) [4]. In this case, the device 
occupies a finite space, the interior of the unit sphere, 



4 



and the refractive index ranges from 1 at the mirror to 2 
in the center. The same trick has been applied in non- 
Euclidean cloaking [l6j |. The mirror creates the illusion 
that light propagates beyond the unit sphere, whereas in 
reality it is reflected at the mirror. After another reflec- 
tion, the light focuses at the mirror image of the original 
focusing point ([3]). To show this, we essay employing the 
inversion in the unit sphere as a mirror transformation 
of the spectator points, 




or, in spherical coordinates, r — ► r . The mirror-image 
of the image ([3]) would appear at 

r-o = -ro ■ (27) 

The inversion (|2"r?)) also preserves the line element (TIT)) 
of Maxwell's fish eye and hence the transformed scalar 
Green's function remains a solution of the wave equation 
(HHJ). We only need to transform the spectator indices of 
the Green tensor by the matrix P (but not the source 
indices). As reflection gives rise to a phase shift of tt [l7| 
we substract the reflected field from the original one, 

G' = G(r)-PG(r- 1 ). (28) 

In Cartesian coordinates, 

P = r 2 l-2r®r (29) 



but in spherical coordinates 

P*,=diag(-r 2 ,l,l) . (30) 

Consequently, at the unit sphere the tangential compo- 
nents of the Green's function, representing the electric 
field, vanish, which is the defining property of a perfect 
mirror. The Green's function (|28|) obeys not only the 
wave equation (J5|) but also the correct boundary condi- 
tions at the spherical mirror, which justifies formula l|28[). 

Finally, we may also include absorption in our theory, 
although in a simplified model. Absorption is described 
by the imaginary part n" of the refractive index. Suppose 
n"(f) is proportional to the real part given by the fish- 
eye profile @. This case is equivalent to having a com- 
plex wavenumber k in the definition (O of the Green's 
function. The singularities of the Green's function (|2"3")) 
describe source and image, but they are not affected by 
the wavenumber that only reduces the amplitude: the 
imaging quality is resistant to absorption, at least in our 
simple model. The fish-eye mirror in 3D images with a 
resolution no longer limited by the wave nature of light. 
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